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J> . Abstract 

in ■ 

\f} \ The scaling dimensions of large operators in M = 4 supersymmetric Yang-Mills 

theory are dual to energies of semiclassical strings in AdSs x S 5 . At one loop, the 
dimensions of large operators can be computed with the help of Bethe ansatz and can 
be directly compared to the string energies. We study finite-size corrections for Bethe 
states which should describe quantum corrections to energies of extended semiclassical 
strings. 

■ 1 Introduction 

J^j ■ The semiclassical regime of the AdS / CFT correspondence |2] relates classical string theory 
^ ■ in AdSs x S 5 to perturbative J\f = 4 supersymmetric Yang-Mills (SYM) theory. This obser- 
vation is quite remarkable and is at first glance counterintuitive, since the "Planck constant" 
of the AdS string theory 1/vA is large at weak 't Hooft coupling A = gy M N in the SYM. The 
key point is that the true parameter of the semiclassical expansion for a large class of string 
states is not 1 / y/X but 1/L [3 , where L is a global charge of the string or some other large 
quantum number. At a given order in 1/L, 1/y/X combines with L into the BMN coupling 
X/L 2 P] that can be large or small depending on L. It is customary now to identify the limit 
of small BMN coupling with the perturbative regime in the dual SYM theory. 

The simplest quantities that can be directly compared to string theory are scaling di- 
mensions of local operators. They are dual to energies of the string sates in AdSs x S 5 . The 
loop corrections to scaling dimensions can be readily computed by standard diagrammatic 
techniques for simple operators, but to describe semiclassical string states we need operators 
with large quantum numbers. Straightforward diagrammatic calculations for such opera- 
tors are hopeless because of the operator mixing which becomes more and more involved 
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as the size of the operators grows. Fortunately, the problem simplifies considerably in the 
large- N limit when only planar diagrams contribute and can be reformulated in a useful way 
which makes it tractable. This reformulation also provides a nice physical interpretation of 
the operator mixing by identifying the mixing matrix with a Hamiltonian of a certain one- 
dimensional spin system, which, quite remarkably, turns out to be integrable and exactly 
solvable [U 03 [7| . The eigenvalues (anomalous dimensions) can then be calculated using 
the Bethe ansatz [HJ |§] and compared to string energies. This has led to spectacular tests of 
the AdS / CFT correspondence for a variety of string solutions with large quantum numbers 



The relationship between classical strings and operators can be established at one and 
two loops quite generally, at the level of effective actions jTHl I2H I2H HZ] or by identifying 
integrable structures j2H|, which are present both in the spin chain jU] and on the string 
world-sheet |2EJ EH EDI- The purpose of the present paper is to go beyond the leading 
order in 1/L in the SYM, which on the string side should correspond to taking into account 
quantum corrections. 1/L corrections to the point-like (BMN) string states pQ were studied 
in E21 122] • They can be compared to anomalous dimensions of near-BPS operators known 
exactly in the first few orders of perturbation theory j^j. The situation is more complicated 
for the extended string solutions. 1/L correction to the string energies are known for a 
particular class of solutions but the dual anomalous dimensions have been so far 

calculated only in the strict L = oo limit. 

We shall consider the subsector of operators of bare dimension L made of two complex 
scalars: Z = $1 + and W = $ 3 + i$A- 



Linear combinations of such operators are in one-to-one correspondence with the states in 
the Hilbert space of a periodic spin-1/2 chain, where Z is identified with spin up and W 
with spin down on each site of the lattice: 



The chain is periodic because of the cyclicity of the trace, which also imposes the condition 
of translation invariance (zero-momentum condition) on admissible states. The one-loop 
mixing matrix in this sector has a form of the Heisenberg Hamiltonian jj]: 



The higher-loop generalizations of this Hamiltonian, potentially also integrable, have been 
discussed in the literature I2H 1213 HHH EE]- The one- loop Hamiltonian can be explicitly 
diagonalized by Bethe ansatz jHJ |§] ■ The Bethe ansatz works also at two and three loops 
[20J and possibly extends to higher loop orders [2E] , but we will concentrate on the one-loop 
anomalous dimensions in this paper*. 

*The Bethe ansatz of (201 ESI is only asymptotic and requires L to be large, but corrections to it seem to 
be exponentially small in the large L limit. 
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The eigenstates of the Heisenberg Hamiltonian with M down spins and with L and M 
large are dual to strings localized at the center of AdS and rotating on the five-sphere with 
two angular momenta J\ = L — M and J 2 = M. If the momenta satisfy 

miJ 1 +m 2 J 2 = 0, (1.2) 

where mi and m 2 are integers, the classical solutions are very simple [SZj ■ The angles on S 5 
are then linear functions of the world-sheet coordinates. We shall concentrate on the gauge 
theory duals of these simplest uniform solutions. 



2 Bethe ansatz 

The vacuum of the Hamiltonian (jl.lj) is an empty state with all spins up and corresponds to 
the chiral primary operator tr Z L . The operators of the form tx{Z L ~ M W M + permutations) 
correspond to the states with M flipped spins and are characterizes by rapidities Ui of M 
magnons. The Bethe ansatz imposes a set of algebraic equations on the rapidities [HIE]: 

Uj+l/2 \ L _yr U j -U k + i 

The cyclicity of the trace requires that the total momentum is zero or, in terms of rapidities, 
that 

nUn + i/2 . 
ipift = < 2 - 2 > 

3 J 



The anomalous dimension is given by 

X 

hi 2 ^ u 2 + 1/4 ' 



-=^£^777- (2.3) 



These equations completely determine the spectrum of the Hamiltonian (jl.lj) . 

We are interested in the limit when L — > 00 and M — > 00 and the filling fraction a = Mj L 
is kept finite. Some particular solutions of Bethe equations in this scaling limit have been 
discussed in condensed matter literature [3B1 EH] and it is this limit that corresponds to the 
AdS duals of the semiclassical strings ^3]. An inspection of the Bethe equations shows that 
Bethe roots scale with L as Uj ~ L. This allows us to simplify Bethe equations a bit. Writing 
them in the logarithmic form, 



Uj + i/2^ f Uj — Uk + i 

and expanding in 1/ui, we get 



Lin " J ^ V = Y In ^±^1 _ 2™, (2.4) 



— + 2t™ = 2^ . (2.5) 
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An arbitrary phase 2irn in (j2.4j) reflects the multivaluededness of the logarithm. In principle, 
different phases can be assigned to different Bethe roots, but a macroscopic number of roots 
should have the same phase to get a meaningful scaling limit at L — > oo. The simplest 
case of equal phases corresponds to the SYM duals of the uniform string solutions [23]. We 
can indeed easily check that (jl.2|) is satisfied. Adding together the rescaled Bethe equations 
()2.5|) . we find the total momentum: 



M 

P 



1 u 
i=i 



but ()2.2|) implies that the total momentum is an integer multiple of 2tt, so 

Mm 

a= — = - (2.7) 
L n 

which is equivalent to (jl.2J) if we put mi = —m, m 2 = n — m. 

Since Bethe equations are even in Ui, the corrections to ()2.5|) are of order 1/L 2 and thus 
()2.5|) is accurate up to and including 0{1/L). In the strict L = oo limit, ()2.5|) becomes a 
singular integral equation, which can be solved explicitly by the Riemann-Hilbert method [23J 
in the general case when the phases are different for different Bethe roots. If all phases are 
equal, the Riemann-Hilbert problem has an algebraic solution. As discussed in the appendix, 
the Riemann-Hilbert approach can then be generalized to take into account the leading 1/L 
correction. Here we use another method based on the observation of J3U] (essentially present 
in the earlier work |U H21 SHI HI] ) that the exact solution of ()2.5|) can be expressed in terms 
of the roots of associated Laguerre polynomials. 

Let us briefly review the arguments of jlU]- Consider the function defined as 

M 

Q{u) = \[{u-u k ) (2.8) 

k=l 

where {wj} are the roots of (|2.5|) . This function is known as the eigenvalue of the Baxter 
Q-operator [IS]. The derivatives of Q evaluated at a root u i: satisfy 



UiQ"(ui) - (L + 2TrnUi)Q'(ui) = 

in virtue of ()2.5|) . and thus the function uQ" \u) — (L+2nnu)Q' \u) is a polynomial of degree M 
which has the same roots {u^} as Q{u). Hence, this function is just Q(u), up to a coefficient. 
Comparing the w M -terms one finds that 

uQ"{u) - (L + 2nnu)Q'(u) + 2nnMQ(u) = 0. (2.9) 

The polynomial solution of this differential equation is the associated Laguerre polynomial: 
Q(u) cx L A j L+1 \2'Knu). So, the roots of the set of equations ()2.5|) are the roots of 

L M {L+I) (27mwi) = 0. (2.10) 
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Some precautions have to be made here since the upper index of the associated Laguerre 
polynomial is negative. Using the Rodrigues representation and the fact that L + 1 > M we 
find: 

L-r^) - ^ (ff - d-«" ( m~-\ ) §• 

We can now use this result to calculate the anomalous dimension in ()2.3|) using the 
relation 

yl = ±I^9M (212) 

where T is a curve encircling all {ui} counterclockwise. We replaced uf + 1/4 by it 2 , because 
Ui = O(L) and we are interested only in the 1/L correction to the anomalous dimension. 
The only singularity outside of the contour of integration is a pole at u = . The contour 
can be deformed such that the integral picks up the residue: 

T\ = - res (W>) = -(2«nf ( L <°' L "'°> - { L '^) . (2,3) 
^fu} ^-»WQ^)J \ (L(0)) 2 ) 

Using the explicit form of the Laguerre polynomials (|2.11j) . we find: 

L(0)L"(0) - (L'(0)) 2 M(M-l) M 2 a{a - 1) 



(L(0)) 2 " L(L-l) L 2 L — l 

where the definition of the filling fraction (|2.7|) was used. Finally, (|2.3|) yields: 



(2.14) 



7 . ^ (l + i j + ( y _ ^ (l + i ) + ( h ) ^ ( , 15) 

This is our final result. The finite-size correction to the anomalous dimension turns out to 
be surprisingly simple and is just proportional to the leading order. 



3 Discussion 



We computed the finite-size correction to the energies of the simplest Bethe states that 
describe anomalous dimensions of large operators in M = 4 SYM. It would be interesting 
to calculate finite-size corrections to anomalous dimensions for a general Bethe state. In- 
deed, the scaling solution of the Bethe equations is known in complete generality and can 
be expressed in terms of Abelian differentials on hyperelliptic Riemann surfaces [2H]. The 
resolvent of Bethe roots plays the central role in this construction, and one may hope that 
calculations in the appendix can be reformulated in a similar geometric language and gen- 
eralized to include arbitrary Bethe states. The calculations in the main text are based on 
another method, which reveals an interesting connection to the Baxter's Q-operator. The 
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eigenvalues of the Q-operator are known to obey a finite difference equation equivalent to 
the set of Bethe equations [IS]. On the other hand, the same eigenvalues satisfy a linear 
differential equation ()2.9|) in the scaling limit. We were not able to derive (j2.9|) directly from 
the Baxter equation, but it is possible that some relationship between them exists. 

The string duals of Bethe states discussed in this paper are known at the classical level 
and are described by very simple string configurations in AdSs x S 5 . It would be very 
interesting to compute quantum corrections to their energies and to compare them to the 
1/L corrections to anomalous dimensions calculated in this paper. However, the relevant 
string solutions are unstable p77j . and it is not quite clear how to take into account quantum 
corrections for them''". This problem is rather puzzling, since the Bethe states at hand seem 
to be perfectly well-defined eigenstates of the Heisenberg Hamiltonian and show no signs of 
instability. The resolution of this problem may provide further insights on the relationship 
between TV = 4 SYM and strings. 
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Appendix A: Finite-size corrections from loop equa- 
tion 

The scaling limit of the Bethe equations (j2.5J) has the same form as the saddle point equation 
for eigenvalues of a Hermitian random matrix j^B]- The analogy becomes complete if we 
rescale Bethe roots and define Xi = Ui/L, which is finite in the limit of L — ► oo. The rescaled 
variables satisfy 



An efficient way to solve matrix models, which turns out to be useful also in the present 
context, is to reformulate the problem in terms of the resolvent 




(A.l) 



and 




(A.2) 




M 



(A.3) 



i'We are grateful to A. Tseytlin for a discussion of this point. 
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The resolvent can be regarded as a generating function for eigenvalues of the commuting 
Hamiltonians of the Heisenberg model [T^ I17j: 

P = -G(0), 7 = -^ Z G"(0), ... (A.4) 

In the matrix models the resolvent satisfies the loop equation ^3 BE] > an analog of which 
can be derived by multiplying both sides of (jA.l|) by l/(x — x^) and summing over %. Then 

2 ST 1 1^1/1 1 

/ ^ Z—^^ [ rp . nn . \ I /y ry> .1 / £1 ' J nn . nn . \ nn nn . nn ry . 

i—t , _ \ ? J / \ % ) ' ■ i ■ % J \ ^ J 

1 v 1 - 1 V 1 1 V 1 

L 2 (x — — a^-) L 2 ^r' — Xi)(x — xj) L 2 [x — Xi) 2 

= G 2 (x) + yG'(x). 

Similar manipulations with the left-hand side give: 

xG 2 (x) + yxG'(x) = G(x) + 2irnxG(x) - 2-Km, (A.5) 

where the momentum condition (jA.2j) in the form (jA.4|) was taken into account. We can 
now expand the resolvent in the powers of 1/L: 

G{x) = G {x) + jG^x) + O (J^j (A.6) 

and plug this expansion in (JA.5j) : 

L° : = xGl - (1 + 2vmx)Go + 2vrm (A.7) 
L- 1 : = 2xG G 1 - (1 + 2nnx)G 1 + xG' . (A.8) 

The leading order is an algebraic equation whose solution is 

G = — ( 1 + 2nnx - J (I + 2nnx) 2 - Mnax] , (A.9) 
2x V / 

where f|2.7|) has been used to express m in terms of the filling fraction a. Alternatively, we 
could derive (|2.7|) from the loop equation by imposing the boundary condition G(x) — > a/x 
at infinity. The first correction can be found from equation ([A.8|i : 



G — X ^'° 



1 + 2nnx — 2xGq 



1 

2x 



1 + 27rn(l - 2a)x 1 



1 + 2imx) 2 - 8imax ^/(l + 2nnx) 2 - 8nnax 
Thus, we have the resolvent up to order L -1 . Its Taylor expansion 



(A.10) 



G(x) = 2-nna + (2nn) 2 a(a -l)(l + -)ac + ... \. I I 



and (|A.4|) now yield 

Xn 2 a(l -a) / l\ 

in agreement with (|2.15|) . 
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